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We calculate the branching ratios and CP asymmetries of the B — > <f)K decays using 
perturbative QCD approach, which includes fcr and threshold resummations. Our results of 
branching ratios are consistent with the experimental data and larger than those obtained from 
the naive factorization assumption and the QCD-improved factorization approach. 



I. INTRODUCTION 

CO 

£f) \ Recently the branching ratios of the B — > cj)K decays have been measured by the BaBar M, BELLE M 

^*>0 and CLEO collaborations. These decay modes are important to understand penguin dynamics, because the 

r B — > (pK modes are pure penguin processes. 

There is a following problem related a penguin contribution for amplitude |4| . A naive estimate of the loop 
diagram leads to P/T ~ a s /(12Tr) \og(m 2 /m 2 c ) ~ 0(0.01) where P is a penguin amplitude and T is a tree 
amplitude. But experimental data for Br(_B — » Ktt) and Br(_B — > 7T7t) leads to P/T ~ O(0. 1). The penguin 
contribution is enhanced dynamically. This enhancement can not be explained by factorization assumption for 
nonleptonic two-body B meson decays B . This problem is explained by Keum, Li and Sanda using perturbative 
^; QCD (PQCD) approach @. 

-— PQCD method for inclusive decays was developed by many researcher over many years, and this formalism 

J> , has been successful. Recently, PQCD has been applied to the exclusive B meson decays, B — > Ktt ||, 7T7t [Q, Tip, 
ttuj M, KK M and Krf ) |lffl. We calculate the branching ratios for the B — > (f>K modes using this approach to 






x 

leading order in a s and leading power in 1/Mb where Mb is the B meson mass. The B — » (f>K modes depend 
only on penguin amplitudes, and the branching ratio for this modes is proportional to \P\ 2 . Thus the B — > tfiK 
modes are useful to seeing the size of the penguin amplitudes. 

In this paper we compute the branching ratio for the B — + tf>K modes using PQCD approach. In Sec. ||, we 



introduce the PQCD formalism for exclusive B meson decays. In Sec. Ill, we calculate analytic formula of the 
factorizable amplitudes and the nonfactorizable amplitudes of the decays. In Sec. [V, we predict the branching 
ratios of the decays and the direct CP asymmetry of the charged mode. Our predicted branching ratios 
agree with the current experimental data and are larger than the values obtained from the naive factorization 
assumption and the QCD-improved factorization mW fll2| ] . In the Appendix, we review the meson wave functions 
with up to twist-3 terms which are proportional to moi^/Ms or M^/Mb- 
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II. PQCD FACTORIZATION THEOREM 



The effective Hamiltonian for AS = 1 transitions is 
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where V* and V q b are the Cabibbo-Kobayashi-Maskawa matrix elements |13|] and Oi_io are local four-fermi 
operators. The local operators are given by 
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where «, j is the color indices and q is taken u, d, s and c. The Wilson coefficients 0i_io are calculated up to 
next-to-leading order by Buras et al. (wl . 



PQCD approach is based on the three scale factorization theorem [[L5| |16| in which the matrix element can 
be written as the convolution of the hard part H with the hard scale t ~ O(Mb), the meson wave function 
$ with the soft scale 1/6 ~ Aqcd and the Wilson coefficient C whose evolves from the W boson mass Mw 
down to the scale t. b is the conjugate variable of the parton transverse momenta kr- For instance, the B — > K 
transition form factor is written as 
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(3) 



where x is the momentum fraction of parton and 7^ is the anomalous dimension of meson. The hard part H 
can be calculated perturbatively and the Feynman diagrams are described as Fig. [j] in leading order of a s . 
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FIG. 1. The B — > K transition form factor in leading order of a s 



It is important that the scale of the Wilson coefficient is determined dynamically, since the scale depends 
on x and b. Then in PQCD, the scale of the Wilson coefficient can reach lower scale than M^/2 which are 
usually taken in the naive factorization assumption case. The Wilson coefficients for the penguin operator are 
enhanced as the scale evolves to lower scale. Therefore, the penguin amplitudes are larger than those in the 
naive factorization assumption. 

In this method, we can calculate not only factorizable amplitudes but also nonfactorizable and annihilation 
amplitudes which can not be calculated in the factorization assumption. 



This B — > K transition form factor has two types of double logarithmic corrections in infrared divergences. 
The one is a s In (Pb) from the overlap of collinear and soft divergences which are generated by the correction 
of the meson wave function as Fig. g(a) and Fig. g(b). The resummation of this double logarithms leads to 
the factor exp[— s(P, b)] JL7| p8[ . The explicit form of this factor is shown, for example, in Ref. [E9[. This 
resummation is called fcy resummation. The other double logarithm is a s ln (1/x) from the end-point region 
of the momentum fraction x [|20|| . This double logarithm is generated by the corrections of the hard part as 
Fig. ||(c) and Fig. ||(d) . The resummation of this double logarithms is called the threshold resummation. We 
use the approximate form of this resummation factor that is T(x) — N t {x(l — x)} c where N t = 1.775 and 
c = 0.3 pij] . The Sudakov factor from kx and threshold resummations suppress large b region and end-point of 
x where the hard part H is singular, the hard part H can then be calculated perturbatively. 
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FIG. 2. These are the diagrams which generate double logarithms corrections in infrared divergences. Fig. H(a) and 
Fig. 0(b) are the correction of the meson wave function. Fig. 0(c) and Fig. g(c) are the correction of the hard part. 

In PQCD approach, we can calculate hard part perturbatively, and we use the meson wave functions which 
are formed in the light-cone QCD sum rules. The hard part depends on processes, but the wave functions are 
independent of those. We can determine the parameters in the wave functions, since there are many modes in 
the B meson two-body decays. So, we can predict quantitative values for the other decays. 

In this study, we calculate the hard part in leading order of a s . The meson wave functions are used with 
leading twist and twist-3 terms which are proportional to ttiok/Mb or M^/Mb- We review the meson wave 
functions in the Appendix. 



III. B -k(>K AMPLITUDES 

In the light-cone coordinate, the B meson momentum P±, the K meson momentum P^ and the 4> meson 
momentum P3 are taken to be 

p 1 = ^|(i,i,o T ), p 2 = ^|(i-^ 0;Ot); p 3 = ^| (r 2 ;lj0T); (4) 

where r^ = M^/Mb, and the K meson mass is neglected. In this situation, the B meson is at rest. The quark 
and anti-quark momenta in the mesons are displayed in Fig. pi Their components are 

fci = (O.ajiPf.kir), Pi-ki = {Pf, (1 - sci)Pf ,-ki r ) , 

fc 2 = (x 2 P+,0,k 2T ), P 2 - k 2 = ((1 - a; 2 )P 2 + ,0,-k 2T ) , 

fc 3 = (0,a;3Pr, k 3T), P3 - k 3 = (0, (1 - x 3 )P^, -k 3T ) . (5) 

We can integrate over kf and choose kf = 0, since the hard part is independent of k± . 
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FIG. 3. Momentum assignment 



The <j> meson longitudinal polarization vector e^ and two transverse polarization vector c^t are given by 

1 
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V2r, 



-(-rll,0 T ), e 0T = (O,O,l T ). 
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The B meson wave function for incoming state and the K and <f> meson wave functions for outgoing state 
with up to twist-3 terms are written as 
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[M ^4> 4 ,{x 3 )+ U /P3^(x 3 ) + M^(x 3 )] , 



(7) 
(8) 
(9) 



where i and j is the color indices, a and /3 are the Dirac indices, took is related to chiral symmetry breaking 



scale, m K = M\j{rad + m s ). v and n are dehned as v^ 



P$/P£ and n» 



z %/ z - = (0,1, T ). We neglect 



the wave functions which are proportional to the transverse polarization vector eT, because these terms do not 
appear in our calculations. It must be noted that the expression of $L°" depends on how to define the sign of 
e^. The explicit form of these wave functions will be shown in Sec. |V] and the detail is shown in Appendix. 
The decay rates for the B — > <f>K have the expressions 
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The decay amplitudes A and .A corresponding to B° — » (^.FC and B° — > 0^° respectively, are written as 

a = UV ts v t tF? + v ts v; b M P + f B v ts v; b F p + v ts v; b M p , (11) 

A = UV t * s V tb F e p + V*V tb M P + f B V t * s V tb F p + V*V tb M p . (12) 

F e is the amplitude for the factorizable diagrams which are considered in the factorization assumption. F a and 
M. are the annihilation factorizable and the nonfactorizable diagrams which are neglected in the factorization 
assumption. The indices e and a denote the tree topology as Fig. ||(a) and annihilation topology as Fig. ||(b) 
respectively. The index P denotes the contribution from the diagram with a penguin operator. The decay 
amplitudes A + and A" corresponding to B + — > <f>K + and B~ — > <pK ~ respectively, are written as 



a + = UV ts v* b F p + v ts v; b M P + f B v ts v t * b F p + v ts v; b M p - f B v us v: b Fj 
a- = f4>v; s v tb F e p + v; s v tb M P + f B v; s v tb F p + v*v tb M p - f B v; s v ub F^ 



v us v: b M T a , 
v: s v ub M T a , 



(13) 
(14) 



where there are very small contribution from tree diagrams, the index T denotes these contribution. 



A. Factorizable Amplitudes 
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FIG. 4. Feynman diagrams contributing to factorizable amplitudes 
The factorizable amplitudes which come from Fig. are written as 

F e = F 3e + F ±e + F 5e 

/>1 />oo 
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x | [x^iix^ixs) - 2r K r^ p K (x 2 ) ((<^(z 3 ) - <t>%{x 3 )) - x 3 (<^,(x 3 ) + <^(x 3 )))] 

xE ai {t^)N t {x 3 {\ ~ x 3 )YK{l - x 2 ,x 3 , b 2 , b 3 ) 

- [(1 - x 2 )<f>£(x 2 )M x 3) + 2r K r <t> (2^(x 2 ) - x 2 0£(x 2 ) - (f>%(x 2 ))) ^(0:3)] 

x J B a4 (i^)iV t {x 2 (l-x 2 )} c / lQ (x3,l-a; 2 ,63,6 2 )| , (16) 

F a p 6 = 16ttC f M| / cfo 2 cfo3 / b 2 db 2 b 3 db 3 
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X< [2^^(x 2 )^(x 3 ) - r^x^ix-i) ((f>l{x 3 ) - (t>%{x 3 ))] 

xE a6 {t^)N t {x 3 {l - x 3 )YK{\ - x 2: x 3 , b 2 , b 3 ) 

+ [r K (l ~ x 2 ) {4>k(x2) + 4> T K {x 2 )) 4>^{x 3 ) + 2r^{x 2 )4>%{x 3 )] 

xE a6 (t^)N t {x 2 (l - x 2 )Yh a (x 3 ,l - x 2 ,b 3 ,b 2 )\ . (17) 



The tree amplitude Fj is the same as F p 4 , but the Wilson coefficient is different. N t {x(l — x)} c is the factor 
for the threshold resummation. N t is a normalization factor and c is a constant, we use N t = 1.775 and c = 0.3. 
The evolution factors are defined by 



E e {t) = a s (t)a e (t) exp[-S B (t) - S K (t)} , E ai (t) = a s (t)a ai (t) exp[-S K (t) - S^{t)} . 



(18) 



The explicit forms of the factor 5,(4) are given, for example, in Ref. g. The hard scales t, which are the typical 
scales in hard process, are given by 



4 1} = max(V^M B , I/61, 1/62) , 4 2) = max( A /xIM s , l/b lt l/b 2 ) , 
i« = max(V^M B , I/62, 1/63) , ti 2) = max( % /T _ ~x^M B , l/b 2 , 1/63) 



(19) 



The Wilson coefficients are given by 
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where a J is the Wilson coefficient for the amplitude Fj . The hard functions, which are the Fourier transfor- 
mation of the virtual quark propagator and the hard gluon propagator, are given by 
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B. Nonfactorizable Amplitudes 
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FIG. 5. Feynman diagrams contributing to nonfactorizable amplitudes 
The nonfactorizable amplitudes which come from Fig. |5| are written as 



M? = M^ 4 + M^ + M P 
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X< [(1 -Xl -X 3 )4>j c {x 2 )-r K X 2 {(/>k(X2) - 4>k{X2))] KSd )h d (xi,x 2 , x 3 ,h, 61,63) 
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p 3 is added to M. P A term, because Aif 3 is the same as Aif 4 without the Wilson coefficient. 
M P =M P 4 + M P 6 
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The tree amplitude A4„ is the same as 
defined by 

KM = a s (t)a> et (t)cxp{-S(t)\ b2=bl } , E' al {t) = a s (t)a' al (t)cxp[-S(t)\ b2=bi 
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The hard scales t are defined by 
Ai) _ 



f = max(£>M B , J\D?\M B , l/b x , l/6 3 ) , tj? = max(FM B , yJ\F?\M B , 1/h, I/63) 
The Wilson coefficients are given by 
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The hard functions are given by 
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IV. NUMERICAL RESULTS 



We summarize the parameters which is used in the numerical analysis of our calculation. 
We use the model of the B meson wave function written as 



4> B (x,b) = N B x 2 (\ - x) 2 exp 



f fxM B 
2 I w B 



where lu b = 0.40 GeV p2| , N B is determined from normalization condition Eq. (B5) 
The K meson wave functions are given as 
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where px = (Md + m s )/Mj( J23| [p24| , The parameters of these wave functions are given as a\ = 0.17, a-i 
0.20, 773 = 0.015, uj 3 = —3.0 where renormalization scale is 1 GeV. We fix the parameters as above values. 
The 4> meson wave functions are given as 
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We fix the parameters as above values, since the numerical results is 



where (3 = 0.024, 6+ = 0.46 |_ 
insensitive to these parameters. 

We use the Wolfenstein parameters for the CKM matrix elements A — 0.819, A = 0.2196, i?t = ^/ p 2 + rj 2 
0.38 [Em, and choose the angle (f>3 = n/2 roj. In addition, we use the following parameters. 



M B = 5.28 GeV, M K = 0.49 GeV, M4 
/ B = 190MeV, /x = 160MeV, h 



= 1.02 GeV, m 0K = 1-70 GeV, 
= 237 MeV, /J = 215 MeV, 



t b o = 1.55 x 10 12 sec, t b ± = 1.65 x 10 12 sec 



A^ D = 0.250 GeV. 



We show the numerical results of each amplitude for the B° — » (f>K° and B^ — > (pK^ decays in Table I. 
The factorizable penguin amplitude F^ is dominant contribution to the B — ► (j)K decays. The factorizable 
annihilation penguin amplitude F^ generates large strong phase. In the B — > 4>K ± modes, there are the 
contribution from JbFJ and M^. These tree amplitudes contribute a few percent to the whole amplitude, 
since the CKM matrix elements related to the tree amplitudes are very small. 

We predict the branching ratios for the B — > <f)K decays as 
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The current experimental values are summarized in Table O. The values which are predicted in PQCD are 
consistent with the current experimental data. However, our branching ratios have the theoretical error from 
the 0(a 2 s ) corrections, the higher twist corrections, and the error of input parameters. Large uncertainties 
come from the meson decay constants, the shape parameter wb, and ttiqk- These parameters are fixed from 
the other modcs(S — ► Ktt, Dtt tttt, etc.). We try to vary lob from 0.36 GeV to 0.44 GeV, then we obtain 
Br(B ± -» 0if ± ) = (7.54 - 13.9) x 10~ 6 . Next, we set uj b = 0.40 and try to vary m 0K from 1.40 GeV to 1.80 
GeV, then we obtain Br(S ± -> 0if ± ) = (6.65 - 11.4) x 10" 6 . 

In the naive factorization assumption, the branching ratio is very sensitive to the effective number of colors 
N° s . If we set N° =3, then the branching ratio is about 4.5 x 10~ 6 where the scale of the Wilson coefficient 
is taken to Mb/2 and F BK is 0.38 from the BSW model. Our predicted values are larger than those obtains 
from the naive factorization assumption. This is due to the enhancement of the Wilson coefficient. In PQCD 
approach, the scale of the Wilson coefficients, which is equal to the hard scale t, can reach lower values than 
Mb/2 in case of the naive factorization assumption, and the Wilson coefficients are then enhanced [||. Therefore, 
our branching ratios are enhanced over those obtained from the naive factorization assumption. 



In the QCD-improved factorization, the branching ratios for the B — > cf>K decays are predicted as Br(S° — * 
(f>K°) = (4.0±?:|) x 1CT 6 and Br(B~ -+ #f~) = (4.3±?;2) X 10" 6 with annihilation effect @. Our predicted 
values are larger than these values. 

We also compute the direct CP asymmetry defined as 



,,lir \A | — \A 

\A+\ 2 + \A~v 



A CP - I A + \2 , I /I -12 ' ( 45 ) 

If the amplitude of the B ± — » 0if ± decay can be written as 



4+ = 1/^e^p - V as V; b e i5T A T , (46) 

A- = V*V tb e i5p A P - V: 8 V ub e i5T A T , (47) 

where the indices T and P denote the tree and penguin respectively, 5t and Sp are the strong phase, and the 
amplitude At and Ap are real, then the direct CP asymmetry is given by 

2X 2 R b sin 03 sm(6p - 5 T )^- 
A%p = \ J ■ ( 48 ) 

Our results of the direct CP asymmetry is about few percent. The reason why the direct CP asymmetry is very 
small is that the CKM matrix element for the tree amplitude is much smaller than for the penguin amplitude. 



V. SUMMARY 

In this paper, we calculate the B° — > <pK° and B^ — > tpK^ decays in PQCD approach. Our predicted 
branching ratios agree with the current experimental data and are larger than the values obtained by the naive 
factorization assumption and the QCD-improved factorization. This is because the Wilson coefficients related 
to penguin operators are enhanced dynamically in PQCD. 

In the PQCD approach, we can calculate the strong phase from nonfactorizable amplitudes and annihilation 
amplitudes explicitly. We can then calculate the direct CP asymmetry of this modes, however this asymmetry is 
very small since the CKM matrix element for the tree amplitude is much smaller than for the penguin amplitude. 

Note added: 

After this work has been completed, we become aware of a similar calculation by Chen, Keum and Li |27|| . 
Our results are in agreement. 
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APPENDIX A: WILSON COEFFICIENTS 

The Wilson coefficients for the m b = 4.8 GeV are 
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Ci(m 6 ) = -2.693 x 1CT 1 , C 2 (m b ) = 1.123 , 

C 3 {m b ) = 1.243 x 1(T 2 , C 4 (m fc ) = -2.665 x 1(T 2 , 

C b {m b ) = 7.753 x 1(T 3 , C 6 (m b ) = -3.256 x 1(T 2 , 

C 7 {m b ) = 3.426 x 1CT 4 , C s (m b ) = 3.131 x 10~ 4 , 

C 9 (m b ) = -9.732 x 1(T 3 , C w (m b ) = 2.217 x 1(T 3 . (Al) 



APPENDIX B: B MESON WAVE FUNCTIONS 

We consider the B meson wave function <&BaB i? = (0\bpj (0)d a i(z)\B (Pi)) where the indices a and (3 are the 
spin indices, i and j are the color indices. This wave function can be decomposed into some components with 
different spin structure. The B meson wave function is then decomposed to be 

*%l M = ^=Jd Xl d*k 1T e-^ p ^ + -^ T ) [( A75 )^( Xl , klT ) +l5<j) P( Xl MT)] a0 , (Bl) 

where N c is color's degree of freedom. The axial- vector wave function ^(aJi,kir) is defined by 

Pi lt Jdxi<Pk.i T e- i ^ p i' + - 1 '^(l4(xi ) ki T ) = 5 -^(O|6(O) 757M d(0)|6(P 1 )) , (B2) 

and the pseudo-scalar wave function <j)^{x\, 1s.it) is defined by 

dsi^kxre-^^i'^-^-'^Csi.kir) = -^=<0|6(0)7s6(«)|d(Pi)) . (B3) 

Using HQET, we can denote <^(xi, kir) = Mb<Pb(xi, 1s.it)- The B meson wave function is then written by 

4*2^ = ■^Jdx 1 <Pk 1T e- i ^ p i* + - k ^ [(A + M B ) 75 B (^,k lr )]^ . (B4) 

The normalization of this wave function is given by 

/ da;i ^(x 1 ,6 1 =0) = -^=, (B5) 

where &i is the conjugate space of ki, and Jb is the decay constant of the B meson. 

APPENDIX C: K MESON WAVE FUNCTIONS 

The K meson wave function is defined that <&K a a ~ (®\s Pj (fydaiiz^K (P 2 )) , z — (0, z~ , Ot) is the coordi- 
nate of the d quark. $jj™^ ■ ■ can be decomposed into some components with different spin structure 757^, 75, 
750>i/, • • •• These matrix elements are calculated using the twist expansion in the light-cone QCD sum rules. 
The twist expansion is the expansion in terms of mojc where m^K = Mj^/(md + m s ). The leading twist terms 
are independent on tuqkj the twist-3 and twist-4 terms are proportional to ttiqk and rn^ K respectively. Up to 
twist-3, the matrix elements are given as [E3| [E4j 
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(0\s(0) l5 ^d{z)\K(P 2 )) = -t/xft M I dxie 

Jo 



■ix 2 P 2 -z 



(j>v{xi) 



(0\s(0)l5d(z)\K(P 2 )) = -if K m 0K / dx 2 e- lx ^ z cj )p (x 2 ) , 

Jo 

/ i\ t 2 \ r 1 

|s(0)75^ d (^)l^( p 2)) = J/ K moif ( 1 - —f- ) (iV* - JW„) / cfe 



-ix 2 P 2 -z 



4>a{X2) 



1 f d 

^K'm,QK{v il n v -v v n v ) I dx 2 e~ lX2P2 ' z - — </> CT (x 2 ) , 
6 p J eb 2 



(CI) 
(C2) 

(C3) 



where v^ = P 2 / ' P 2 and n M = z^/z = (0,1, Ot)- 4>v{x 2 ) is the leading twist wave function, (f) p {x 2 ) and 
d/dx 2 (f) a (x 2 ) are the twist-3 wave functions. These wave functions are normalized as 



/ dx 2 cj) v (x 2 ) = 1,1 dx 2 cj) p (x 2 ) = 1 , / dx 2 - — <j> a (x 2 ) = . 
Jo Jo Jo dx 2 



We define the following new wave functions, 



<t>K^ 2 ) 



/. 



K 



--4> V (X 2 ) , <1>k( x 2) 



I 



K 



2V2N C v ' " ' 2V2iV ( 
The K meson wave functions are then written as 



4p{%z) , 4>k{x 2 ) 



I 



K 



12j2N~ c dx 2 



(j)a{x 2 ) 



$ 



(in) _ —iS- 

K.afi.ij 



*7 



f 2N c J 



-ix 2 P 2 -z 



(C4) 



(C5) 



[/fWs^O^) +™ x75^(a;2) + m 0K (/h /O - l)j 5 <f>K(x 2 )] a . (C6) 



APPENDIX D: <£ MESON WAVE FUNCTIONS 



The (j> meson wave function has much terms than the B meson and K meson wave functions, since the <f> 
meson is a vector meson. The wave function $\ ™K .. = (0\s/3j(0)s a i(z)\<j)(P3)} can be decomposed into some 
components with different spin structure 7^, 7^75, & au , I, ■ ■ ■. The components are calculated using the twist 
expansion in the light-cone QCD sum rules. In this case, the twist expansion is the expansion in terms of M^, 
Up to twist-3, the components are given as p5| 



(0|*(0)7^«(«)|^(ft)) = UM* / dx 3 e~ lx ^ z 



^^*//(*s) + e^(*,) 



= UAlJ dx 3 e- ix > p °- z ^u<j>//(x 3 ) + eT^\x 3 ) , 
Jo L J 

[0\- S (0) 7 , l5S (zM(P 3 )) = -\ ^U ff-^j M^eT v P 3a z p J\x 3 e^ p ^g { °;\x 3 ) 

= -\UM^; a ^l v P i<x Z J dx 3 e'^ p ^g^(x 3 ) , 
(O|s(O)^ S (z)|0(P 3 )) = ifj f dx 3 e- ix ° p °- z \ (e^P 3 „ - e^P^) c/> ± (x 3 ) 



+ (P3uZ„ - P^zJjg^M$hf)(xs) 



Jo 



(P3-Z) 2 
{t^Psv ~ ^ P 3ij) &±{x 3 ) 

(e4.fj.P3v - e<t>vP3ii)h/j (X3) 



(Dl) 



(D2) 



(D3) 
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= -fjM^J 1 dx 3 e~^ p ^l^-h ( ;/(x 3 ) , (D4) 

where n^ = z^/z + ,v^ = P 3ti /P 3 ~. <p//(x 3 ) and (j>±(x 3 ) are the leading twist wave functions, h, (x 3 ), -j^-hu (x 3 ), 
9± (^3) an d <?_l (#3) are the twist-3 wave functions. These wave functions satisfy the normalization conditions 



I dx 3 (p//(x 3 ) = 1 , / dx 3 hj/(x 3 ) = 1 , / t£r 3 - — /*y(a:3) = 0, 

1 ,.1 ,.1 

da; 3 </ij.(x3) = 1 , / efo30X (^3) = 1 , / dx 3 gY'(x 3 ) = 1 . (D5) 

Jo Jo 

We can neglect the wave functions which are proportional to the transverse polarization vector eT, because 
these terms do not appear in our calculations. We define the following new wave functions, 

M*s) = t4t^3) , #(*») = JLrrfjte) > n^) = X^//fe) ■ (D6) 



The meson wave functions are then written as 

^p,tj = -^= J dx 3 e-**>*-' [M UM^)+ A &<%(xa) + M^%{x 3 )] . (D7) 



It must be noted that this expression is dependent on how to define the sign of e^. 

We considered the wave functions for an initial state. In this study, we need the K meson and <p meson wave 
functions for a final state. The wave function for a final state is given as $( ou *) = 7o$^™^7o- 
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B° -> fiK 

-1.03 x 10" 1 
6.45 x 10~ 3 + i 4.28 x 1CT 2 
5.24 x 10~ 3 - i 3.61 x 10~ 3 

-8.03 x 10~ 4 - i 1.73 x 10~ 3 



B i 



■(j>K d 



UF e p 

f B F p 

M P 

M p a 

f B F T 

MT 



-1.03 x 10 _1 

6.17 x 10~ 3 + i4.20 x 10~ 2 

5.24 x 10 -3 - i 3.61 x 10~ 3 

-6.56 x 10~ 4 - i 7.22 x 10~ 4 

-1.11 x 10 _1 - i 3.75 x 10~ 2 

1.60 x 10~ 2 + i 2.77 x 10~ 2 



TABLE I. Contribution to the B° — > cf>K° and B — > <j>K decays from each amplitude. 



Br(B° -» <f>K°) 

5.1+25 ±0.8) x 10 ~ 6 

!.7J;f;g ± 1.5) x 10" 6 

< 12.3 x 10" 6 



Bt(B ± -* cf)K d 



BaBar 

BELLE 
CLEO 



(7.7tiA ± 0.8) x 10" 6 
(10.6+ 2 ;g±2.2) x 10~ 6 
(5.5J: 2 ;8 ± 0.6) x 10~ 6 



TABLE II. The experimental data of the B -> <j>K branching ratios from the BaBar jjj], the BELLE Q and the CLEO 
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